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1 Introduction 



We investigate the following Cauchy problem: 

{dtu{t,x)+d:,{i^){t,x)+l[u{t,.)]{x)-dl,u{t,x)=Q iG(0,r),xGM, 
[ti(0,x) = uq{x) x G M, 

where T is any given positive time, uq G L^(M) and X is a non-local operator defined as follows: for any 
Schwartz function Lp G and any x G M, 

A^]{x):= / |c|-3^"(x-CK. 

Remark 1. Equation ([T]l can also be written in conservative form 

dtu + ^ Y + C[u] - dxu^ = 

where 

m{x)--= / \(:\-^^'{x-c)dc 



Equation ([TJ appears in the work of Fowler mm on the evolution of dunes; the term dunes refers 
to instabilities in landforms, which occur through the interaction of a turbulent flow with an erodible 
substrate. Equation ([B is valid for a river flow (from left to the right) over a erodible bottom u{t, x) with 
slow variation. For more details on the physical background, we refer the reader to ||5ll3. 



1 



Roughly speaking, I[u\ is a weighted mean of second derivatives of u with the bad sign; hence, this 
term has a deregularizing effect and the main consequence is probably the fact that O does not satisfy 
the maximum principle (see below for more details). Nevertheless, one can see that the diffusive operator 
—d^^ controls the instabilities produced by Z and ensures the existence and the uniqueness of a smooth 
solution for positive times. The starting point to establish these facts is the derivation of a new formula 
for the operator I, namely This result allows first to show easily that Z — d^^ is a pseudo-differential 
operator with symbol ipxiO — ^vr^^^ — aj|^|3 + i ^xCl^l^' where aj and hj are positive constants 
(see (111)). The symbol Air"^^^ corresponds to the diffusive operator —d'^^ and — aj|^|3 + i fej^lCl^ is 
the symbol of the nonlocal operator Z. Notice that this last symbol contains a fractional anti-diffusion 
— axICI^ (recall that this is the symbol of —{—d^^)'^, up to a positive multiplicative constant) and a 
fractional drift i 6j^|.^|3. Because of the fact that the fractional anti-diffusion is of order |, the real 
part of il^xiO behaves as up to a positive multiplicative constant, as ^ ^ +oo. A consequence 
is that Equation ([T]) has a regularizing effect on the initial data: even if uq is only L^, the solution u 
becomes for positive times. The uniqueness of a L°^((0, T); L^) solution is obtained by the use of 
a mild formulation (see Definition [T]) based on Duhamel's formula (fT2)) . in which appears the kernel K 
of Z — d'^^. The use of such a formula also allows to prove local-in-time existence with the help of a 
contracting fixed point theorem. Such an approach is quite classical; we refer the reader, for instance, to 
the book of Pazy IH and the references therein on the application of the theory of semigroups of linear 
operators to partial differential equations. We also refer the reader to the work of Droniou et al. in 131 
for fractal conservation laws of the form 

dtu + d^{f{u)) + {-dl,)^[u] = Q, (2) 

where / is locally Lipschitz continuous and A G (1, 2], and to the work of Tadmor ||9l on the Kuramoto- 
Sivashinsky equation: 

dtu + ]^\d.M^-dl,u = {-dl,f[u]- 

In fact, fractal conservation law (jj)) is monotone and the global existence of a L°° solution is based on the 
fact that the norm of u does not increase. In our case, this is not true and we have to use a classical 
energy estimate to get a global estimate. The regularizing effect on the initial data are first proved by 
a fixed point theorem on the Duhamel's formula to get regularity in space and next by a bootstrap 
method to get further regularity. This technique has already been used in |j3J. 

On the other hand, one of our main result is probably the proof of the failure of the maximum 
principle for ([T|l: more precisely, we exhibit positive dunes which take negative values in finite time, 
since we establish that the bottom is eroded downstream from the dune. We also give some numerical 
results that illustrate this fact (for more precision, see Remark|2]and Section|7]). The proof of the failure 
of the maximum principle is based on the integral formula Q. Roughly speaking, this formula means 
that T is a Levy operator with a bad sign, see [2 |. Notice that the Kuramoto-Sivashinsky equation is also 
non-monotone, but no proof of the failure of the maximum principle is given in [9|. 

The paper is organized as follows. In Section |2j we give the integral and pseudo-differential formula 
for Z; we also establish the properties on the kernel K ofZ — d"^^ that will be needed. In Section [3l we 
define the notion of mild solution for ([T]). Sections |4] and [5] are, respectively, devoted to the proof of the 
uniqueness and the existence of a mild solution; Section|5]also contains the proof of the regularity of the 
solution. The proof of the failure of the maximum principle is given in Section [6l Finally, we give in 
Section |7] some numerical simulations that illustrate the theory of the preceding sections. 

Here are our main results. 
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Theorem 1. Let T > and uq G L^(M). There exists a unique mild solution u G L°°((0, T); L^(M)) o/ 
(O f Definition [7]). Moreover, 

i) u £ C°°((0, T] X M) and for all to G (0, T], n and all its derivatives belong to C([to, T]; L^(IR)). 

//j n satisfies dtu + dx{^) +I[u] — d'^^u = 0, o?i (0, T] x M, /n classical sense (I[u] being 
properly defined by ^ and ©j. 

Hi) u G C([0, T]; L^(]R)) and n(0, .) = almost everywhere (a.e. for short). 

Proposition 1 (L^-stability). Let {u,v) be solutions to ([T]l with respective initial data {uo,vo), we 
have: 

\W - ^llc([0,T];L2{R)) < C {T,M,\\uo\\l2(k.),\\vo\\l2(r)) \ \uo - Vo||l2(R) 

where M := max (| |u| |c;([o,t];L2(M)), I l^^l Ic([o,t];L2{M)))- 

Theorem 2 (Failure of the maximum principle). Assume that uq G (M) n (M) is nonnegative and 
such that there exist G M with uo{x^) = Uq{x^) = Uq{x^) = and 

Then, there exists t* > with M(t*, x*) < 0. 

Remark 2. Hypothesis of the theorem above are satisfied, for instance, for non-negative uq G C^(M) H 
//^(M) ™c/i f/zaf exists £R. with uo{x^) = u'q{x^:) = Uq{x^) = and 

Vx < X*, uo(x) > and 3xo < x* s.t. no(xo) > 0. 

A simple example of such an initial dune is shown in Figure [7] Observe that the bottom is eroded 
downstream from the dune (recall that the nonlinear convective term propagates a positive dune from the 
left to the right). 



Figure 1: Evolution of a dune, at i = and t = t^,. We can observe that u{t^:,x^:) < and that 
/ u{t, x) dx remains constant. 



Notations: In the following, we let denote the Fourier transform defined for / G L^(M) by: for 
all ^ G M, 

•^/(O := / e-''-^^fix)dx. 



We also let T define the extension of the preceding operator from to . In the sequel, we only 
consider Fourier transform with respect to (w.r.t. for short) the space variable; in order to simplify the 
presentation, for any u G C([0, T]; L2(M)), we let Tu G C([0, T]; L2(M, C)) denote the function 

t G [o,r] G L2(M,C). 



2 Preliminaries 

In Subsection 12.11 we give the integral and the pseudo-differential formula for X and in Subsection 
we give the properties on the kernel of J — 

2.1 Integral formula for J 

Proposition 2. For all Lp G 5(]R) and all x G M, 



-oo 



,|7/3 



with Ci = |. 



Proof. The proof is an easy consequence of Taylor-Poisson's formula and Fubini's Theorem; notice that 
the regularity of tp ensures the validity of the computations that follow. We have: 

/■O ip{x + z)-^{x)-^'{x)z /-O . .7 f ^ „, ,2^\j 

J — ^Jy^3 ^ ' dz = J \z\ l^J {l-T)ip"{x + Tz)z^dTjdz, 

= J {1-t)(^J \z\~^ip"{x + Tz)dz^ dr, 



j\l - t)t-1 (^IJ^ \C\-h"ix - C)dc) dr, 



thanks to the change of variable tz = Then, 

ip(x + z) — (p(x) — u:>'(x)z , , , 2 . , , , , 

^ |^|7/3 dz = (1 - r)T-sdTl[ip]ix) = -iMx). 

The proof is now complete. ■ 
Corollary 1. There are positive constants aj and bj such that for all ip G S{R) and all ^ G M, 

:F{i[^]-^"){0 = MO^m, (4) 

where ipiiC) = 47r^^^ - ax\^\t + i 
Proof. We have 



oo 



1 7/3 
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Notice that Proposition |2] ensures that for Lp G G L^(M) and thus its Fourier transform is 

well-defined. By Fubini's theorem, we can first integrate w.r.t x to deduce that 

where we let T^z^ denote the (translated) function x — > ip{x + z). Classical formulae on Fourier 
transform imply that (^) = ■ip{£,)Tip{^), where 

Simple computations show that 

It is immediate that the real part of ip{^) is even, non-positive, non-identically equal to and homoge- 
neous of degree | (the last property can be seen by changing the variable by z' = ^z). Moreover, the 
imaginary part of -0(0 is odd, 
constants aj and bx such that 



3 

imaginary part of -0(0 is odd, negative and homogeneous of degree | on . There then exist positive 



and, in particular-, (0 = (-axlCl^ + i hm^' ) ^^iO- Since = ^tt^C^J^^{0^ 

the proof of Corollary [T]is complete. 
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Remarks. 1. Since I[ip] = 1r^\ ■ {'s * ip" ,we have J='{I[ip]) = J='{1rJ ■ {'s) ■ {-4tt^\^\^) ■ J=-{(p). 
Elementary computations give .F(1]r_^| • |~3) = r(|) — isign(^)-^^ |'^|~^- Hence ax = 
T{l)\andhj = T{l)^. 

2. Let s(zM.Ifip£ H^{M.), one can also define 2[(p] through its Fourier transform by 

Hmm) ■■= -4vr2r(^) -^sign(0^^ 1^1^ -^M 

Thus, if if E H^, we have that I[ip] G H^^i and \\T[ip\\\ ^^_a < A7r'^T{^)\\ip\\H^ ■ This implies in 
particular that I : H'^{R) C;,(M)nL2(M), since by Sobolev embedding ^ Cb(M)nL2(R). 

3. Corollary\!\implies that I - dl^ : C'^{W) n if2(M) ^ C(]R) n L'^{M) with 1 which satisfies both 
formula ^ and Q. 



2.2 Main properties on the kernel K oil — dl^ 

By Corollary [H we see that the semi-group generated by T — d"^^ is formally given by the convolution 
with the kernel (defined for i > and x G M) 

K{t,x)=J'~^{e-''^^){x). 
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Proposition 3. K{t, •) is a real valued continuous function. 

Proof. K{t, •) is a real valued continuous function as inverse Fourier transform of a W'^'^ function 
with an even real part and an odd imaginary part. ■ 

In the sequel, we only consider real valued solution of ([IJ. We expose in Figure |2] the evolution of 
K{t, •) for different times. Note that K{t, •) is not compactly supported but that K{t,x) < for 

\x\>iv,ithc{t) = ^m.nK{t,.)m\\L^ 




Figure 2: The kernel of T - for t = 0.05, 0.1 and 0.5 s. 



Proposition 4. The kernel K has a non-zero negative part. 
Proof. Let us assume that K is nonnegative, then 

= I K{t, .)=J' h'-^ie^''^^)) (0) = e-*^^(°) = 1 

4 § 

for all C G IK; hence, since je-^'^^^^)] = e-*^^'''!^!'-"^!^!^) > 1 for < |C| < this gives us a 
contradiction. ■ 

The main consequence of this is the failure of the maximum principle for the equation 

dtu+l[u]- dl^u = d] (5) 

that is to say, there exists a non-negative initial condition uq such that, for some t > 0, u{t, .) := 
K(t, .) * uq has a non-zero negative part, see section |6] below. Nevertheless, K enjoys many properties 
similar than those one satisfied by the kernel of the heat equation and that ensure that Equation ^ has 
a regularizing effect on the initial condition: if uq G LP(R) for some p e [1, +oo), then u is C°° for 
positive times, see section [5] 

Let us precise here the properties that will be needed in this paper. Since K(t, •) G (M), the family 
of bounded linear operators {uq G L^(M) — > K{t, •) * no G L^(M)}^^q is well-defined. Moreover, it is 
a strongly continuous semi-group of convolution, that is to say: 

Vt, s > 0, K{s, •) * K{t, •) = K{s + t,-), . 
Vno e L2(M), limt^o K{t, ■)*uo = uq in L'^{R). ^ ' 
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Next, the kernel K is smooth on (0, +00) x M and we have: 

vr>o, 3/Cos.t. vt G (o,r], ||a,i^(t,.)||L2(u) </Cot-^, (V) 

VT>0, 3/Cis.t. Vt G (0,T], ||a,i^(t,.)||Li(M) </Cit-i (8) 
Vt, s > 0, K{s, •) * d^K{t, •) = d^K{s + t, •). (9) 

Proof of these properties. The semi-group property Q and Q are immediate consequences of Fourier 
formula. Let us prove the strong continuity. By Plancherel's formula, 

\\K{t,-)*uo-uo\\l2^R-j = \ \J='{K{t,-)*uo) - J='uo\\l2(^^) 

= We-'^^J'uo - ^no||i2(K) = / \e-'^^ - 1|' I-FuqI'. (10) 

The function |e~*'^^ — Ip |.FnoP converges pointwise to on R, as t — > 0. Recalling that minRe{ipx) is 
finite, we infer that le'^'l'^ - Ip | J"uoP < C| J"moP and the dominated convergence theorem implies that 
the last term of (ITOl) tends to as t ^ 0. This completes the proof of (O. Let us now prove the estimates 
on the gradient. The smoothness of K is an immediate consequence of the theorem of derivation under 
the integral sign applied to the definition of K by Fourier transform. We get in particular: 

Since the function 2m^e 

-ti>i{0 is d^K{t, •) is L2 and we have: 

Jr Jr 
Let us change the variable by ^' = 1 2^. We get: 

Jr 
Jr 

for all t G (0, T]. The proof of Q is now complete. To prove ([U, we have to derive a "homogeneity-like" 
property for K. Easy computations show that 

Jr 

^ f g2i7rx5g-t(47r2|5|2-ai|5|i+i6ie|5|^)^^^ 

Jr 

by changing the variable by ^' = 1 2^. Then, 

^(^^3.) ^ ^-1 /" g2i^(t-^x)5'g-(4^2|?'|2-ai|5n3+ifei«'IS'l^)e-(l-i^)(ai|5'|3-ibi€'|r|3)^^/^ 
JR 

= /" e2i7r(f-^x)|'g-,/.l(C')g"(l--t3)(ai|S'|3-i6i^'|5'|^)^^/_ 
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Fort < 1, define G((l - := jr-i |^e^(i-t^)(ai|5'|^-i . It is readily seen that G is 

as inverse Fourier transform of a W^'^ function. Moreover, for Iq G (0, 1) and all t G (0, to]> 



\mi-ts),-)\\LHR)<c 



< C{to), 



where C(to) only depends on to- Classical formula on Fourier transform then give: 

K{t, x) = t-5 [k{1, ■) * G((l - rs), .)) {t--2x). 

Observe now that d:,K{l, •) = ^ 2i ^vre-^^^^)) is as inverse Fourier transform of a ' 

function. Then, 

d,Kit,x)=t-^(^d,Kilr)*Gi{l-d),-)) 
is L and its L norm can be computed by the change of variable x = t 2 x as follows: 

||9,ir(t,0llL^(M) =t-^P.K(lr)*G((l-ti),OllLi(R) <*-^P.K(lr)llLi(M^ 
for any t G (0, to]- Since 

P,K(t,-)||ii(M) <C||e^2ievre-*'^^(«)||^2.i(R,c) <C(to,r), 
for all t G [to , , the proof of ^ is now complete. ■ 
Remark 4. For anj no G (M) anJ t > 0, 

||K(t,-)*no||i2(R) <e""*||no||L2(u), (11) 

where cuq = — minRe(?/'i). 

Proof. This is readily established with Plancherel's formula, like in (fTOl) . ■ 

3 Duhamel's formula 

Using Fourier transform and Corollary [H we formally see that any solution to ([B satisfies Duhamel's 
formula (fT2l) (see also the proof of Lemma[3l which justifies the computations). This observation is the 
starting point of the definition of mild solution below. 

Definition 1. LetT > and uq G L'^{R). We say that u G L°°((0, T); L2(M)) is a mild solution to ([B 
iffora.e. t G (0, T), 

1 /■* 

u{tr) = K{t,-)*uo- - d^K{t- s,-)*u^{s,-)ds. (12) 
^ Jo 

The following proposition shows that all the terms in ([T2l ) are well-defined and that Equation ^ 
generates a (non-linear) semi-group. 
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Proposition 5. Let T > 0, uq £ L'^{R) and v G L°°((0, T); L^(M)). Then, the function 

1 /■* 

u:te{0,T]^K{t,-)*uo-- da:K{t-s,-)*v{s,-)dseL'^{R), (13) 

2 Jo 

well-defined and belongs to C([0, T]; L^(]R)) (being extended at t = by the value u{0, .) = uq). 
(Semi-group property) Moreover, for all to € (0, T) all t & [0,T — to], 

1 /■* 

u{to + *,.)= •) * ^(*o, •) - o / "^x-fs^l* - s, •) * ^(*o + s, ■)ds. 

^ Jo 

Proof. By it is classical that the function t € (0, T] —>■ K(t, ■) * uq £ L^(]R) is continuous and can 
be continuously extended by the value u{0, ■) = uq att = 0. What is left to prove is thus the continuity 
of the function ^ 

w:te[0,T]^ / d^K{t- s,-)*v{s,-)ds e L^iR). 
Jo 

Let us extend dxK and v for all times the following way: 



n{t, 

Then we have 



if not ^'^10 if not. 



w{t,-) = / n{t- sr)*V{s,-)ds. 
Jr 

It is immediate that V G L°^(M; L^(IR)). Moreover, Q imphes that 



and it follows that TC £ L^(M; ^^(R)). Young's Inequahties imply that for all t £ 



(14) 



\\n{t-s,-)*V{s,-)\\L2(^K)ds < / \\n{t-s,-)\\L2^u)\\V{s,-)\\LHmds, 

Jr 

< \\'^\\l'^{R]L2{r)) \\^\\l°°{R;L'^{r:))- (15) 
This implies, in particulars that the function w is well-defined. Let us now take t, s G M and define 



/ nit-T,-)*V{T)dT- [ n{s-T,-)*V{T)dT 

Jr Jr 



We have 



/ < / ||(7^(^-r,•)-W(s-r,•))*V(r)||^.(JJ)dr, 

< / ||7i(i-r,.)-W(s-r,.)||L2(i,)||V(r)||ii(M)dT, 
Jr 

thanks to Young's Inequalities. It follows that 

I< \\n{t -T,-) -n{s -T,-)\\L2(^)dT\\V\\Lo^(^.Li(Ryy 

Jr 
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Since the translation are continuous in L^(]R; L^(M)), we see that / ^ as |t — s| — > 0. In particular, 
the function w is continuous and this completes the proof of the continuity of u. 
Let us now prove the semi-group property. By ^ and (O, we infer that 

1 Z"*" 

u{to + t,-) = K{t,-) * K{to,-) *uo - - dxK{t + to - s,-) *v{s,-)ds 

^ Jo 



t+to 

dxK{t + to - s, •) * v{s, ■)ds, 

to 

to 



K{to, •) * K{t, •) * no - J / ° K{t, •) * d^^Kito - s, •) * v{s, ■)ds 

^ Jo 

1 /■* 

-- d^K{t-s',-)*v{to + s',-)ds', 



thanks to the change of variable s' = s — tQto compute the last integral term. Then, 



to 



uito + t,-) = K{t,-)*K{to,-)*uo- K{t,-)*^ I d^K{to-s,-)*v{s,-)ds 



- 1" 



dxK{t — s , ■) * v{to + s' , ■)ds' , 

dxK{to - s,-) *v{s,-)ds^ 

d^K{t - s', •) * v{to + s', ■)ds', 

i f d.,K{t-s\-)*v{to + s\-)ds'. 
^ Jo 



K{t,-)* [K{to,-)*uo 2 



= K{t,-)*u{to,-) 2 

The proof of the semi group property is now complete. ■ 
Remark 5. For v € L°°{{0,T); L^{R)), u € C{[0, T];L'^{R)) defined in ^ satisfies: 

IfIIc{[0,T];L2{R)) < e'^" IFo||l2(k) + 2/C0T4 ||ti||^oo((o,T);Ll(M))- (16) 

Proof. Indeed, with ([T4l ) and ([T5] ). we estimate the integral term of (fT3] ) and with ([TTI ). we estimate the 
norm of K{t, ■) * uq. ■ 



4 Uniqueness of a solution 



Let us state a lemma that will be needed later. 



Lemma 1. Let T > 0, uq £ L'^{R). For i = 1,2, let Vi G L°°((0, T); L^(]R)) and define Ui £ 
C([0, T]; L2(]R)) as in Proposition^by 



1 



Ui{t, •) := K{t, ■)*uo- - / dxK(t - s,-)* Vi{s, ■)ds. 

^ Jo 

Then we have the estimate 

Ikl - 'U2||c([0,T];L2(lK)) < 2K,qT~^ \\v\ - 1 |l°°((0,T);L1 (M)) ■ (l^) 
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Proof. For all t G [0, T], we have 



1 r 

ui{t, ■) - U2{t, ■) = -- / dxK{t - s, •) * (■"1(5, •) - V2{s, ■))ds. 
^ Jo 



Hence, 



\ui{t, •) - U2{t, 



< 



dxK{t - s, •) * {vi{s, •) - V2{s, ■))ds 



\d^K{t - s,.)* {vi{s, .) - V2is, .))IIl2(r) ds. (18) 



By ©, 

\\da,K{t -s,.)* {vi{s,-) - W2(s,-))IIl2(r) < \\dxK{t - s,-) ||L2(iR)||t;i(s, •) - V2{s, OIIlHk) 

3 

< ICo{t - s)~i\\vi{s,-) - V2is,-)\\LH^y 

Inequality (fTSl ) then gives 

KL 3 

•) -'"2(i, •)||l2(R) < J^it- s)~^ds\\vi-V2\\L'^((^o,ty,Lim), 

= 21Cot^ \\vi - V2\\L^(^(o,t)■,L^{R))■ 
ln paiticular, for all s E [0,t\ 

\\ui{s,-) -'U2(s,-)||l2(ir>) < 2/CoS3 11^1 - t'2||L-'((0,s);Ll(M)) < 2/Cot4 \\vi - ^^2 1 ((0,t);Ll (R)) 

and we have proved that 

- 1 1 

\\U1 - U2\\c{[0,t];L^R)) < S/Cot* \\vi - ^2 1 1 L°° ((0,t);Ll (R)) • (19) 



Proposition 6. Let T > and uq G L^(M). There exists at most one u e L°°((0, T); L^(M)) w/z/c/i a 
mild solution to O. 

Proo/ Letu,v G L°°((0, T); L2(]R)) be two mild solutions. Lett G [0,r]. With Lemma [U applied to 
vi = and V2 = v^, we get 



k - ^llc([0,t];L2(E)) < 2/Cot4 \\u^ - U^||Loo((o,t);Li( 



(20) 



Since ||^^^-'y^||Loo((o,t);Li(R)) < M||M-t;||c([o,t];L2(iR)) withM = | |w| |c([o,t],l2{r)) + | bl Ic{[o,t],l2(r)), 
we get: 

Ik - 'f^llc{[0,t];L2{R)) < 2M/Cot4 11''^ - ^^||c([0,t];L2(R))- 

We then have established that tt = u on [0,t] for any t £ {0,T] such that t < (2M/Co)"^. Notice that 
since u and w are continuous with values in L^, u = v on [0, T*] with = (2Af/Co)~'^ > 0. To prove 
that u = f on [0, T], let us define to := sup{t G (0, T] s.t. n = f on [0, t]} and let us assume that to / T. 
The continuity of u and v implies that ti(to, ■) = ^(to, •)• The semi-group property of Proposition |5]thus 
implies that u{to + ■,■) and t;(to + •, •) are mild solutions of ([Hi with the same initial condition; that is to 
say n(to + 0, •) = 'y(to + 0, •). The first step of the proof then implies that n(to + •,•)= f (to + •, •) on 
[0, m.m{T^,T — to}]; hence, we get a contradiction with the definition of to and we deduce that to = T. 
The proof of the uniqueness is now complete. ■ 
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5 Existence of a regular solution 

This section is devoted to the proof of the existence of a solution u £ C^'^((0, T] x ]R)to ([T]); that is to 
say, u is in space and in time. We first need the following technical result: 

Lemma 2. Let uq G L^(M) W T >0. Letve C([0, T]; Li(R)) n C((0, T]; W^i'i(]R)) that satisfies 

sup ti\\d^v{t,-)\\L^-^^ < +00. (21) 
te(o,r] 

Let u G C([0, T]; L2(M)) be thefi^nction defined in (O. Then, u G C((0, T]; ^^^(M)) with 

1 I 1 1 

sup t2||aa;'u(t,-)||L2{R) < ''Ci||no||L2m) + — - T4 sup 1 2 | , .) 1 1 (22) 

te(o,T] ^ te{o,T] 

1 3 1 

where I is a constant equal to Jq (1 — s)~^s~2ds = B{l/2, 1/4), B being the beta function. 
Moreover, let Vi G C([0, T]; ^^(M)) n C((0, T]; ^^^'^(M)) satisfy (EB W Jf'^ie -Uj by ^ (with u and 
V replaced, respectively, by ui and Vi)for i = 1,2. Then, 

sup t2\\da:{ui -U2){t,-)\\L2(^K) < Ti SUp t2\\d^{vi - V2){t, ■)\\l^R)- (23) 

te(o,T] ^ te(o,T] 

Proof. Recall that Proposition [5] ensures that u G C([0, T]; L^(M)). It is easy to check that the distribu- 
tion derivative of u w.r.t. the space variable satisfies: for any t G (0, T], 

1 /■* 

dxu{t, •) = dxK{t, ■)*uo- - / dxK{t - s, •) * dxv{s, ■)ds. 

^ Jo 

Let us verify that all the terms are well-defined in L^. Since dxK{t, •) G L^(M), it is obvious that 
dxK{t, •) * no G L'^{R). Moreover, define 

1 /■* 

w{t,-)-=- / dxK{t- s,-)*dxv{s,-)ds. 
^ Jo 

Young's Inequalities and ^ give 

\\da:K{t - s, •) * dxv{s, •)||l2(R) < \\dxK(t - S, ■)\\L2(R)\\dxV{s, 

= \\dxK{t - s, •)||l2(jj)s"^s^||9^.v(s, OIIl^r), 

< }Co{t-s)-U-^ sup T^\\dxV{T,-)\\Lim. (24) 

Te{o,T] 

Since Jp(t — s)~3s~2 ds < oo, by (|2T]) we deduce that w{t, ■) is well-defined in and thus for all 
t G (0, T], dxu{t, •) G L^(M). Let us now prove that dxU is continuous on (0, T] with values in L^. For 
5 > Oandt G (0,r], define 



ws{t 



1 /■* 

, •) := / -s,-)* {l{s>s}dxv{s, •)) ds. 

^ Jo 
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Since l{s^sydxv{s, •) G L°°{[0, T];L^{M.)), Proposition |5] ensures that ws is continuous on [0, T] with 
values in L^. Moreover, for any to G {0,T],6 < Iq and t G [^O) 



1 



& 



•) - •)||l2(ir) < -j^ j^\\dxK{t- s,-)*dxv{s,-)\\L2(^^)ds, 



- ^ / {t - sr^s'-^ds sup s^\\dxv{s,-)\\L^R) bydal, 
^ JO se(o,T] 



- ^ [ i^O-s) 4S 2(^5 sup S2\\dxV{s,-)\\Li( 

'0 i*e(o,T] 



2 



It follows that 



/'(^ 3 1 1 

sup \\w{t,-) - Ws{t,-)\\LVR) < ^ (to - s)~^S~^ds sup S2\\dxV{s,-)\\Li(R) ^ 0, 
ie[io,T] ^ Jo s6(0,T] 

as 5 — > 0. We deduce that w G C((0, T]; L^(R)) as local uniform limit of continuous functions. More- 
over, 

dxKit, ■)*uo = T-^ (i ^ 2i T^ie-'^^^<i^Tuo{i)\ . 



The dominated convergence theorem immediately implies that for any to > 0, 

2 



This means that t > ^ — > 2i ■n^e~^'^^^^^ Tu^) G L^(]R) is continuous and, since T is an isometry 
of I?, we deduce that t > ^ dxK{t, •) * liQ G L^(M) is continuous. We then have established that 
dxU G C((0, T]; Let us now estimate how the norm of dxU can explode at t = 0. By (|24l) . 



•)IIl2(m) < ^ / (t-s)"4s"2ds sup r2||9j,7;(T, OIIli 

^ JO re(0,Tl 



= * sup T2\\dxV{T,-)\\Li{M), 

^ re(0,T] 

1 3 1 

where / = /g (1 — s')~^ s'~^ds' = i?(l/2, 1/4); notice that the last integral term has been computed 
with the help of the change of variable s' = f ■ Moreover, ^ and Young's Inequalities imply that 

\\dxK(t, •) * Uo\\l2{R) < /Cit~5||Mo||L2(K). 

We deduce that for any t G (0, T], 

1 J'C I 1 1 

Pxn(t, •)||l2(i.) < /Cit~2||no||i2(K) + — sup S2\\dxV{s,-)\\LH^), 

^ se(o,T] 

which imphes immediately (|22)) . 
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Let us now prove (l23l) . For any t G (0, T], 







/"^ 3 1 1 

< ^ / (t-s)"3s"5ds sup -t,2)(s,-)||il( 
/O sG(0,T] 



2 



-— sup s2||5^(t'i - f2)(s, OIIlmiR)' 
^ se(o,T] 



which implies immediately (l23l) . 



Remark 6. Lef uq,T,v and u that satisfy assumptions ofLemma\2\ Then, we have established that for 
any t G (0, T\, 

1 /■* 

•) = <9a;i<:(t, •) * ^0 " 7^ / dxK{t - S, •) * ■)ds. 

^ Jo 

Let us now prove the local-in-time existence of a regular solution. 

Proposition 7. Let uq G L^(]R). There exists > on/j depends on ||'Uo||l2(r) ^mc/j ([D 
admits a (unique) mild solution u G C([0, T*]; L^(]R)) n C((0, X,,]; ff^(M)) on (0,T^:)such that 

sup t^||(9a;u(i, Olliam) < +00 anJ sup t\\dl^u{t, ■)\\L2m < +00. 

te(o,n] te(o,T,] 

Moreover, u belongs to C^'^((0, T,,,] x M) satisfies the PDE in ([U m f/ie classical sense. 

Proof We use a contracting fixed point theorem. Forn G C([0, T*]; L2(IR))n C((0, T*]; define 
the norm 

llklll := ll^^llc([o,T.];L2(R)) + sup 1 2 1 1 a^n(t , •) | |l2 (R) • (25) 

iG{0,T,] 

Define the space 

X := {u£ C{[0,T^];L'^(M.))nC{{0,T^];H^(R)) s.t. n(0,-) = uq and |||n||| < +00} . 

It is readily seen that X is a complete metric space endowed with the distance induced by the norm 1 1 1 • 1 1 1 . 
For u £ X, define the function 



1 /■* 

Qu:te [0, T*] K{t, ■)*uo- - / d^K{t - s, •) * u^{s, ^ds G L^(M). 

2 Jo 



(26) 



By Proposition 121 Gn G C([0, T^^]; ^^(k)) and satisfies eti(0,-) = uq. Define v := n^. We have 
= 2ud^u. Therefore that v G C([0, T*]; Li(M)) n C((0, T^]; ^^^'^(M)) and that ^ holds true. 
By LemmalU we deduce that @u G X. Let us take R > \\uq\\i2(^^ + /Ci| |tio| |l2(jj) and assume that 
|||n||| < R. Since | |l°°((o,t*);L1{M)) = 1 1""! lc'([o t*]-l2(r))5 estimate (fT6l ) of Remark |5] implies that 

l|0^llc([O,T,];L2(R)) < e'^° *||uo||l2(ir) + 2/Cor/ ||u||p([g^^_^].^2(K)), 

< e"o^*||7Xo||L2(K^ +2/Cor^ i?^ (27) 
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Estimate ((22l) of Lemma 111 implies that 

sup t2\\dx{&u{t,-))\\L2(^) < /Ci||uo||L2m) + — - T*^ SUp \\dx{u^){t, (^), 

tG{o,r*] ^ te(o,TH,] 

~ 2 



by Cauchy-Schwarz inequality. Adding this inequality with (|27] ). we get: 

|||Gn||| <e^o^*||no||L2(K)+/Ci||uo||L2(K) + (2 + /)/Corii?2. 
For G (0, T] sufficiently small such that 

e""^*ll^^o||L2{R) +/Ci||no||L2(K) + (2 + /)/Cor)i?2 < (28) 

we deduce that |||0ii||| < R. To sum-up, we have established that for any T* G (0,T] such that (1281 ) 
holds true, G (defined by (l26l )) maps into itself, where Br denotes the ball of X (endowed with 
the III • III norm) centered at the origin and of radius R. Let us now prove that is a contraction. For 
u,v ^ Bfi, Estimate ([TT] ) of Lemma [U implies that 

1 

\\@U - Qv\\c([0,T,];L^R)) < '^RICqT^' \\u - ^1 |c([0,T.];L2(R)) , (29) 

where we again used | \u'^-v^\\c{[o,T,]■,L^R)) < (M \c{Io,t,],l^r)) + M \c{[o,n],L^R)))\\u-v\ \c{[o,n];L^R)) ■ 
Moreover, Estimate (1231 ) of Lemma E] implies that 

1 i 1 

sup t2\\dx{@U - Qv){t,-)\\L2(^-^-^ < JCqIT^^ sup t2\\{udxU - vdxV){t,-)\\Ll(^^y 
tG{0,T.] tG(0,T,] 

Since 

t^\{ud^u - vd^v){t, •)||li(i8) < t^\\dxv{t, ■)\\l^m)\\{u - v){t, 

+t^\u{t, •)llL2(K)||5x.(n - D)(t,-)||L2(iR), 

< infill ||(n- ?;)(t,-)||L2(iK) 

+llklll t'2\\dx{u-v){t,-)\\L2(^^), 

< i?|||u — u|||, 

we get: sup(g(o,T,] ■)IIl2(r) < RfCnlTJ^ 1 1 |n — f 1 1 1 . Adding this inequality with (1291). 

we find that 

IIIGn- et;||| < (4 + /)i?/CoT)|||u-u|||. 

1 

Consequently, for any > sufficiently small such that (l28l) holds true and (4 + I)RK,qT^ < 1, G is 
a contraction from Br into itself. The Banach fixed point theorem then implies that G admits a (unique) 
fixed point u £ C{[0,n]; L^{R)) D C{{0,n]; H^(R)) satisfying supte(o,T.] OIIl^CR) < oo 

which is, of course, a mild solution to ([T]). 

To prove the regularity of u, we have to use again a contracting fixed point theorem. But, 
this is now the gradient of the solution which is searched as a fixed point. Let to £ (0, T,,). For any 
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t G (0, - to], define u{t, •) := u{to + t, •)■ Let Ti G (0, - to]. We still endow C([0, Ti];L^{R)) n 
C((0, r,']; F^(R)) with the norm ||| • ||| defined in ^ with replaced by Ti. Define the complete 
metric space 

X' := {t; G C([0,r:];L2(R)) n C((0, T^]; s.t. z;(0, •) = vo and |||t;||| < +cx)} , 

where := dxu{0, •). For f G X', define the function 



e'v.te [0,Tl] K{t, ■)*vo- [ dxK{t - s, •) * {uv){s, ■)ds G U 

Jo 



(30) 



Arguing as in the first step of the proof, we claim that Proposition |5l Remark [51 Lemmas [T] and |2] imply 
that ©' maps X' into itself with: for any u,v G X', 

\\\e'v\\\ < e'^°^*\\vo\\mR)+iCi\\vo\\mR) + ctI~^ IMl, 
llie'u-G'u'iii < ct!,^\\\v -w\\\, 

for some nonnegative constant C that only depends on /Cq and ||u^||c([to,T,];_ffi(R))- Let us take R' such 
that 

R' > e'^'''^*\\vo\\L2{R) +^l|bo||L2(K). 

If Ti > satisfies 



e"°^* 1 1^0 1 1 L2(M) +-^11^0^2(15) +CTl^i?' < R' and CT^^ < 1, 

then 0' maps Bjii{X') into itself and is a contraction. Let v denote its unique fixed point. Observe 
now that Q'dxU = dxU, thanks to Remark |6] But, similar arguments than these ones used to prove the 
uniqueness of a mild solution in the preceding section allow to show that there exists at most one function 
w G L°°((0,r^');L2(M)) that satisfies e'w = w. It follows that d^u = v ^ X' on (0,r^'); hence, we 
deduce that u G C(to, to + Tl]-H'^{W)). To sum-up, we have proved that for all to G (0, T*], there exists 
Tl G (0, - to] such that u G C((to, + T'^]; H'^ {R)) . This completes the proof of the continuity of 
u on (0, T^] with values in H^. The proof of the C^'^ regularity is postponed to Lemma|3]in the next 
section, where it will be useful for the maximum principle failure. ■ 

We can finally prove the global-in-time existence. 

Propositions. Letuo G L'^{M)andT > 0. There exists a (unique) mild solution u G C([0, T]; L^(M))n 
C((0, T]; ij2(M)) to © such that 

sup t^\\dxu{t,-)\\L2(^) < +00 and sup t| |a^^n(t, •)! |l2(r) < +oo. 
te(o,T] te(o,T] 

Moreover, u belongs to C^'^((0, T] x M) and satisfies the PDE in ([U in the classical sense. 

Proof. We have to derive first a estimate on the local regular solution u constructed in Proposition |7] 
Multiplying ([T|) by ti and integrating w.r.t. the space variable, we get: 

■^7: [ u^dx + 1 {I[u] — d1xu)udx = 0. (31) 
dt 2 jjj 
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Indeed, the following computations show that the nonlinear term equals 0: 

/ dx (—j udx = — I —dxudx = — / u{dxuu)dx = — / udxi — ) dx. 

But, Corollary [U implies that 

[ {I[u] - dlxu)udx = [ J^-^{il)jru)udx = [ Vx|^up(i^= [ Re(V'j)|J^n|2d^, 
JR Jr Jr Jr 

since Jjg(T[u] — dxxu)udx is real. It follows that, 

/ {I[u] — dxxu)udx > minRe('0i) / \!Fu\'^d^, 
Jr Jr 

= minRe('i/'i) / u^dx, 
Jr 

thanks to Plancherel's Equality. Equation (|3T]) then implies that 



I u^dx < (Jo u^dx 



2 

and by Gronwall's Lemma, we deduce that for all t G [0, Xi,] 

\\u{t, -MIl^r) < e'^''*l|wo||L2(R)- 

Define now 

to := sup{t > s.t. there exists a (unique) mild sol. to ([Hi 

on (0, t) that satisfies the regularity of Proposition [8]} 

and let us assume that to < T (recall that Proposition |7] ensures that to > 0). By Proposition |71 there 
exists > such that for any initial data vq that satisfy ||fo||L2(R) < e'^"*"||no||i2(]g), ([T]) admits a 
regular mild solution on (0,r*) with initial datum vq. Hence, if we define vo = u{to — T*/2), then 
([T]) admits a mild solution v that satisfies the regularity of Proposition [8] Using the uniqueness and the 
semi-group property, it is now easy to show that u{to — + t, •) = v{t, •) for all t £ [0, r*/2] and 
that the function u defined by 5 = u on [0,to] and u{to — + t, •) = v{t, •) for t € [T^,/2,T^,] 
is still a mild solution to ([T]) that satisfies the regularity of proposition [H Since the solution u lives on 
[0, to + ^*/2], this gives us a contradiction. We conclude that to >T and this completes the proof of the 
global existence of a regular- solution. ■ 

Remark 7. To sum-up, we have proved Theorem [7] with the C^'^ regularity of u. To obtain further 
regularity, we claim that we can use the same method by arguing by induction. 

Now we prove the L^-stability stated in Proposition [T] 

Proof of PropositionU}- Let {u, v) be solutions to ([B with respective initial data {uo, vq). let T > 
and t G [0, T]. Substracting 



n(t, 



1 /■* 

K{t, ■)*uo- - / dxK{t - s, •) * u^is, ■)ds 
^ Jo 
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and 



t 

vit, .) = K{t, ■)*vo-^ I d^K{t -s,-)* v\s, ■)ds 



v[t,.) =K[t,-)*vo- ^ 

we get 







1 /■* 

u{t, •) - v{t, •) = K{t, •) * (uo -vo)-- / d^K{t -s,-)* •) - v\s, ■))ds. (32) 

^ Jo 

Hence, by (fTTl) of Remark|4]and Young inequality 

\\u{t,-)-v{t,-)\\L^K) < e'^°'^\\uo-Vo\\LHR) + 2 I l-^^i^ (* " S, •) I |l2 (R) I (s, •) " ^'^ (s, •) 1 1 Li (R) (^S- 

Taking M = max (||u||c([o,t];L2(r)), ||v||c([o,t];L2(r))) , we can bound 

\\u{t,-) - v{t,-)\\L2m < e'^°'^\\uo-vo\\LyR) + M \\dxK{t - s,-)\\l2(^^\\u{s,-) - v{s,■)\\L2r^^ds 

Jo 

< e'^°'^||^/o - ^^o||l2{r) + / {t- s)~^\u{s,-) -v{s,-)\\L2^^-)ds, 

thanks to (O. With lemma |4l the proof is finished. 







6 Failure of the maximum principle 

We now investigate the proof of Theorem |2] We first need a regularity result which ensures that if the 
initial data is regular then so is the solution up to the time t = 0. 

Lemma 3. Let uq G H'^{M.) and T > 0. Assume that u is a mild solution to ([T]) that satisfies the 
regularity of Proposition^ Then, u is in fact C([0, T];H'^(R)) D C^'^((0, T] x M) and satisfies the PDF 
in ([T]) in the classical sense. Moreover, if uq G C^(M), then u G C^'^([0, T] x M) and satisfies the PDE 
up to the time t = 0. 

Proof. First, we leave it to the reader to verify that the continuity with values in up to the time t = 
can be proved again by the use of a contracting fixed point theorem. Note that the regularity of uq allows 
to work in a space of continuous functions with values in up to the time t = 0; more precisely, we 
argue as in the proof of Proposition |7J but we can directly use the C([0, T*]; H"^) norm instead of the 
III • III norm defined in (1251 ). Let us now prove that n is a classical solution to ([IJ. Taking the Fourier 
transform w.r.t. the space variable in ([T2l ). we get: for all t G [0, T], 



jr(n(t, •)) = e-^^^Tuo - fin- e-'^*-'^^^ r{u^{s, ■))ds. (33) 

Jo 



Since G C([0, T]; Li(M)), we know that /"(u^) g C([0, T]; Cb(M, C)). For any ^ G M, the function 
t G [0, T] T{v?{t, G C is thus continuous. Define 



uj{t,i) :=- [% Tr^e-^'-'^'f'^^^^T{u\s,-)){Ods. 
Jo 
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Classical results on ODE then imply that w is derivable w.r.t. the time variable with 



dtw{t,0+MCMt,^) = = --^l 5x ( y ) (t,-) ) (0- (34) 



,2 



Let us prove that all these terms are continuous with values in L^. First, u G C([0, T]; iJ^(M)) therefore 

2 



dx{y?) G (^([0, T]; L^(M)) and we deduce that T{dx{^)) is continuous with values in L^. Moreover, 



Equation (1331 ) implies that 

V'l w{t, •) = •)) - e-'^^Tuo) . 

Since u G C([0, T]; L^(R)) and ipj behaves at infinity as | • ipjw is continuous with values in L^. All 
the terms in (|33] ) then are continuous with values in and this implies w G C^([0, T]; L^(M, C)) with 

|(«;(t, •)) + V'l •) = (y ) (t, •) 

Moreover, it is easy to see that t G [0, T] e'^'^'^TuQ G L^(M, C) is with 

From Equation (l33l) . we infer that is on [0, T] with values in with 

|(^(n(t, •))) = -i^iw{t, ■) - i,je-'^^Tuo - T [d^ (y ) (t, •) 

= -^xT{n(i,-))-T(iudxU){t,-)). 
Since is an isometry of L , we deduce that u G C^([0, T]; L2(]R)) and that 

^Kt,.)) = -a.(y)(t,.)-^-^(V'xW,-))), 

= -5. (y ) (t, •) - X[n(t, •)] + 

where we used Corollary [T] to compute the pseudo-differential term. In particular, u satisfies the PDE of 
([T]) in the distribution sense. What is left to prove is the regularity in space of u. Differentiating (ITll) 
two times w.r.t. the space variable, we get: for any t G [0, T], 



K(t, •) * n'o - I dxK(t - s, •) * v{s, ■)ds. (35) 
Jo 



where v = {dxu)^ + u9^^u. By the Sobolev imbedding i7^(M) ^ ^7^(1^), we know that v G 
C([0, T]; Li(M) n L2(]R)). By LemmaS we know that for all y G M, 

|(9^.ir(t-s, •)*t;(s, •)(2;)-9:rK(t-s, •)*v(s, OC?/)! < \\dxK{t-s)\\L2(^^)\\Ti^x-y){v{s, ■))-v{s, OIIlzcr)- 
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By ©, we deduce that for all t G [0, T] and all x, y G M, 



< 



' dxK{t — s, ■) * v{s, ■){x)ds — / dxK{t — s,-) * v{s,-){y)ds 
Jo 

* _3 1 

/Co(i-s) ■i\\'T^x-y){v{s,-)) -vis,-)\\L2(^)ds <4:Ti sup ||T(3,_y)(t;(s, •)) - ^^(s, •)IIl2( 

sG[0,T] 



By LemmalU we deduce that the second term of (1351 ) is continuous w.r.t. the space variable independently 
of the time variable (equicontinuity w.r.t. the time variable). Moreover, we already know that this term is 
continuous on [0, T] with values in (by Proposition |5]l and Lemma|7]implies that it is continuous w.r.t. 
the couple (t, x) on [0, T] x JR. We now leave it to the reader to verify that {t, x) — > K{t, •) * Uq{x) is 
continuous on (0, T] x M when uq G H^{M.) and continuous on [0, T] x M when moreover uq G C^(M). 
The proof of Lemma [3] is complete. ■ 

The proof of Theorem |2]is now an immediate consequence of the integral formula 

Proof of Theorem^ Lemma|3]and Proposition |2]imply that the solution n to © is C^'^ up to the initial 
time t = and that 

/n ^ , f \ ' f \ , n r ^o(^* + uo{x^) - u'o{x^)z /// n n 

ut[U,x^) +uo{x*)Uq{x^) + C2 / —Yj^ dz-UQ{x^) = 0. 

J —00 \z\ 

It follows that 

n,(0,x,) = -Cx r ^^^%±^dz<0. 
There then exists t* > such that M(t*, a;*) < 0. The proof of Theorem [2] is now complete. ■ 



7 Numerical simulations 

The aim of this part is to show some numerical simulations for ([T]l. An explicit discretization gives results 
in line with the theoretical study (see Remai-k|2l). 

We write ([T]) with a viscous coefficient e > as follows: 

dtu + dx(^Y + -^m) - ^^'-^ = 0' 
where for any G 5(M) and x G M, 

/.+00 

Cy]ix):= / \Cr^ip'ix-C)dC. 



The viscous coefficient is taken sufficiently small, in order to magnify the erosive effect of the non-local 
term. The new definition of the non-local term {I[u\ = dxC[u\) follows 13, which interpretes C[u\ as 
a flow. Notice that in HHH, the bottom is, in fact, s{t,x) = u{t,x + q'{l)t), where q is the bedload 
transport of sediments; for the sake of simplicity, we continue to work with u. 

To shed light on the effect of the nonlocal term, we compare the evolution of the solution of (l36l ) 
with the solution of the viscous Burgers equation: 

dtu + dx (^) - edl^u = 0. (37) 
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7.1 Maximum principle for the viscous Burgers equation 



It is well-known that (ITTl) satisfies the maximum principle: for any initial data uq G (M), ess-inf uq < 
u < ess-sup Mo- As a consequence, (l37l) cannot take into account erosion phenomena. To simulate the 
evolution of u, we define a regular discretization of [0, L] with a spatial step Ax such that L = MAx, 
and a discretization of [0, T] with a time step At such that T = NAt. We let Xj, t„ and uf respectively 
denote the point iAx, the time nAt and the computed solution at the point {iiAt, iAx). We use the 
following explicit centered scheme: 



,n+l 



u? + At 



1 

2 



2ur + uZ 



2Ax 



Ax' 



It is well-known that this scheme is stable under the CFL-Peclet condition: 

'Ax Ax^' 



At = mill 



u 



2e 



(38) 



(39) 



To convince the reader, let us simulate the evolution of the well-known following travelling waves of 
(127]) for e = 1: 



u{t,x) := - 



1 — tanh I ^ f ^ ~ 



We expose in Figure [3]both analytic and numerical solutions. We observe an eiTor of the order of 10 ^ 
between these solutions. Let us now take, as an initial dune, the following small regular perturbation on 




Figure 3: Numerical and analytic travelling waves of the viscous Burgers equation. 



the bottom: 



-1 



uo{x) = \ e-(^-T^ if|_i<^<| + i, (40) 
otherwise. 



We describe its evolution in FigurelH The dune propagates, but as mentioned above the erosion phenom- 
ena are not taken into account since u remains positive (because of the maximum principle). 

Remark 8. Equation ([T]) also admits travelling wave solutions, see [TS. 
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Figure 4: Evolution of the solution of dST]) with uq defined in (gOl {L = 30, M = 4001 and e = 0.1). 



7.2 Erosive effect of the nonlocal term 



Let us return to the study of (|36] ). We add the discretization of the non-local operator C to the explicit 
centered scheme (|38] ). It is natural to consider the following discretization: 



+00 n _ „.n 

y i\ z^u I 2A2; 



Instead we will approximate 



1 wtj+i - 



3=0 



This is based on the assumption that for x € [0, L], 

c[u{t,-)]{x)^ r\cr"sdMt,x-odc. (41) 







This fact is not true for general u, but if we assume that the initial profile uq satisfies uo{x) = 0, Vx < 
and semi-discretize in time Equation (l36l ). we get : 



u{t + At, x) = u{t, x) + At y) - d-^C[u{t, .)] + edl^u^ . 

We observe that u{t + At, x) = 0, Vx < and by induction u{tn,x) = Vx < 0, V?i. Now 

/•X 

C[u{tn,.)]= \C\-^d^u{tn,X-C)dC. 

Jo 
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Actually, we take uo G C^(]R) and supp(tfo) CC (0,L) (see FigurelD. Moreover, Lemma |3] suggests 



The initial dune 




(L/2)-1 (L/2) 



(L/2)+1 



Figure 5: The initial dune defined in (l40l) . 



that all the derivatives of u are continuous with values in w.r.t. the time variable up to the time t 
It then is natural to expect that (at least for small times) Equation (|4TI ) is a good approximation. 
We then use the following explicit scheme for 



uf + At 



1 K+iJ 
2 



\2 



2Ax 



2Ax 



+ e- 



u: 



2< + <-i 



Ax' 



As far as the stability condition, one can numerically see that ( [39l ) is still ensuring stability for small 
Ax. The evolution of the initial dune (l40l ) is given in Picture [6l As the solutions of the viscous Burgers 
equation, the dune is propagated downstream but we now observe an erosive process behind the dune: 
the bottom is eroded downstream from the dune, as shown in Remark|2l 

Let us make a final remark. We are aware of that the fact that these numerical simulations are a first 
crude attempt. To tackle rigorously the non local term would need further study, which will be reported 
elsewhere. 



23 



0,0075 
0.D25 




Figure 6: Evolution of an initial dune, by using the non-local model (l36l ) (L = 30, M = 4001 and 
e = 0.1). 

A Some technical lemmas. 

We first recall a generalization of Gronwall's lemma proved e.g. in lH. 

Lemma 4. Let g : [0, T] be a bounded measurable function and suppose that there are positive 

constants C, A and 9 > such that, for all t <T , 

g{t) <A + c[ {t- sf-^g{s) ds. 
Jo 

Then, 

sup g{t) < CtA, 

0<t<T 

where constant Ct does not depend on A. 

Lemma 5. Let f,gG L^(M). Then, f * g G C(M) and for all x,y gR, 

\f*9{x) - f*g{y)\ < WT^x-y)/ - f\\L^R)\\9\\L^R)■ 
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Proof. The result is immediate if / and g are smooth; indeed, 



f{x - z)g{z)dz - I f{y - z)g{z)dz 



\f * 9{x) - f * g{y)\ = 

< I \f{x-z)-f{y-z)g{z)\dz, 
JR 

< ll'^(x-3/)/ - /IIl2(R)I|5||l2(R)- 

The result for general / and g only L^, is then obtained by density. ■ 

Lemma 6. Letu G C{[0,T]; L'^{R)). Then, supjgp^^j \ \Th{u{t, •)) - u{t, 0, as h ^ 0. 

Proof. The function u is uniformly continuous with values in as a continuous function on a compact 
set [0, T]. For any e > 0, there then exist finite a sequence = to < ti < . . . < tN = T such that for 
any t G [0, T], there exists j G {0, . . . , - 1} with 

\\u{t,-) - u{tj,-)\\L2^u) < 

Moreover, 

\\%{u{t, ■)) - U{t, < \ \%{u{t,-)) -%{u{tj,-))\\L\R) 

+ \\Th{u{tj,-)) - 'U(tj,-)I|L2(R) + \\u{tj,-) - w(i,-)llL2(]R). 

Since \\Th{u{t, ■)) - Th{u{tj, ■))\\l2^^) = \ \u{t,-) - ^(tj, •)! |l2{r)> we get: 

\\Th{u{t,-)) - n(t,-)||L2(R) < \\Th{u{tj,-)) - n(tj,-)||L2(]R) + 2\\u{tj,-) - •)||l2(r), 

< \\Th{u{tj,-)) - n(tj,-)||L2(iK) + 2e. 

By the continuity of the translation in L^(M), \ \Th{u{tj, •)) — u{tj, •)||l2(r) — > 0, as /i — > 0. Then, 

limsup||T/j(n(t,-)) - ^^(i, •)! Il2(r) < 2e. 

Taking the infimum w.r.t. e > implies the result. ■ 

Lemma 7. Let u G C([0, T]; such that u is continuous w.r.t. the variable x uniformly in t. Then, 

u G C([0,r] X M). 

Proof. Let {tQ^xo) G [0,T] x M. Let e > 0. By the regularity of u w.r.t. the space variable, we know 
that there exists r/ > such that for any t G [0, T] and all y G [xq — 7], xq + r]], 

\u{to,xo) - u{t,x)\ < \u{to,xo) - u{to,y)\ + \u{to,y) - u{t,y)\ + \u{t,y) - u{t,x)\, 
< e + \u{to,y)-u{t,y)\+e. 

If we integrate w.r.t. y e [xq — rj, xq + rj], then we get: 
2r]\u{to,xo) -u{t,x)\ < 4er/+ / \u{to,y) -u{t,y)\dy < 4eT?(2r?)2 ||u(io, •) - uit, ■)\\L2^^y 

Jxo-r) 

By the continuity of u with values in L^, 

limsup \u{tQ,xo) — u{t,x)\ < 2e. 

Taking the infimum w.r.t. e > completes the proof. ■ 
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